“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1986-06 


Surface reconstruction from planar contours 


Jones, Allan R. 


http://ndl.handle.net/10945/22134 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
F (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist | | Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


i KNOX appointed — and published — scholarly author. 
wh 


http://www.nps.edu/library 






LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


BARE Aarts Oo HaGid Od. UMMM AR Becbikoto dt CANOE 













































































































4 . 
' . F Lat 94 Sot OO ce Fiat eran 
aa . "Ads ior fe it Aa oh 84a phy KB dd Bete GA Arm i de yak cay Ali tp a, Gch, Aaa led 4.5 da @ 264 fe! Ris Ro erate Sarthe 7] . Sear meats 
| rm seats atin atin atiad oleate haar amen teen Reis 
Parone De to mS; a, 8 ‘ ad t ‘ Uy Be Aetimty 4 aravanes 
A ae athe aia Lhd pike yu sete te ds Oot ie SOs d Bei e Basb he eh Ter) yt Mtr dee ay pratt. AE ver 
naa oa eer. ‘ SORE Wane ds ver &. reir DIANE Xo eM bh bs 45842 By my Me 8 AIRED Ra Be Oe ae 
‘ . b . i. o . oy ae Pa yr et é fh heme VL dtd tb 6st ths ra ery wear boas 9 RY Wi Ale GAR A/D Ble Bem ASAD Mc Midednity Neate 8 
: , ae ia A ‘ 4 ds 4 tee RT eS ; Ch Rad Ore B tabs Mahe yt toe. ages eas, pa a pd Ae yy Vingeras Beaks Mr 5 2: ts ds tales 204 ed 
' ' “2 : Wool Bey ho Gb ge ais AA ta ide BF om t- 2 eee he Mae CAL eV Arteta Ab ore 
‘Ve 5 os f 4 : , - uy hh sf bre parkes bag Me Sarda taal eek d Lethe kos BY mee cde be A BEUSRTA A ty au AZing beg 212 CANOE LAD B20) che. So-m PR Bee 
z . ~ ue Spor Par 4 a aa } ; au r* b} ; ras te) _ ine + wats Noes to OB. Lrbes Ae af US arta wire lielies, vend aulpe Nidetate >: Maus obseds s50s 8.88 Fagg 
' ’ ES eee reals *¢ ° 5 : . : ad pie. BP 4% By 
ae na i 7 es oe a yeriee \! alge ’ ee Bei yb shy ee A UR 1 Shen ao tee Pert ia Ay GES Jia st a ayaa he eivoea semen 
fl (ean ; Site ' $ s » \ . gh a) a state $$ 0 be Y kege per Ute 8d-d ob he yh TAD EB, 8. a a, ara 9 fe Amt Wd pomre eY Det Wnt te 
7 caer pr oath as A ' ars ae earn te | by Ab ae a ae “Lathoric Lae te are ard ae a 20 Ot bebe be %) ihe tytn! fat PE SBS 2, AU Sines Mrnite A188) Be Ber de R-Ar Ro s. 21m Al® ante 
| eC ek Tate ine nes anh OO, atta Bake erties iactataatanc ett arta hte neh ate cadet aan atatagmiopieeeattte 
ahs ; ' ¥ . We ae onde ey oo Pe a er ee Sided te ta deNnied Fs bode 84:9 AAS, Lag We Gave OO ae Mi Acaw or Var sa3le 9 meee 8, bso ¥ as © 
yy 5 @ a A an - tet 3 Reid & Te wt $ PRG ay » bys ‘ ib: eit ace teheg PSUR ot! re oy Ra NE ater Yc edd avis che per an eer.’ Sabre PCT 
c 4 r 4a tra » A vig ; Je e 7a y i Mahe wet med ed Ok rite ey ea a aay ab. bee iy NAT ye WC! gt A at 04 Ot &, 0, A Ak Me WBA TA Ne ar mm O05 My Ms 
’ rs ~ LT t Pe wt ‘ . beim 
: } ; Ae ages ie nowt hy es, ; Sey ee? mee BS kee x A a4 ATi ore aus. aii wi) ed vin Reh eriahetritiy srk ena Ar hn iriirereray apices 
' wt a 8% aula . rie anh " aU N et mb be qeaeee 4 “e tems wae 3 jem . Poe ose ay +. te HM 8. Cat de fea ee tint gles atacn: & As Ro Be 84 1s teh BOR, WeBe A. 4440) 
‘ : a ? 
bn . : a ats Pa ae + mpd RO AA As fT apt. # tt og 4 Wert he 8 40M: 4, 05 SoA TOG: 8 ere iy Tbe dtty ee 4 ta 6d 91. ts ordi ne tne BoB: Oe ss BA Bs 
R c ' ’ i eN2 mint, bev A pea pe Oe ee RO Ae ee eg f a $ va. i a hers 0 Sat Oe : 4 mut a Rehr Bx er be % 4 
; Fae AG rie P » eee Reeve et t x a , ror vi AY MODE Le peilegewa, 1, st viata aca atvanatat gh nah a tase ty shtehe nents 
Ye ees 4 iS : $ 426 are Ooh argh Bt af bal 95 fe Pe & Me Meh Me a he Be Pi . Seas day FA pews Pe A @ 
ae sees : payee oy Ws i Rie of Waris sXs a as Mero 9 oa ere» bege a ee ae weet re, uBa ds Seishin = neon nae veh etree mentee 
i] 2 ‘ we ad “Jha te ai Rattle 4 RN 1? ale wA ~ . . rhe +f bens 
: Sie att ZUG TORU Eh Fh ~~ Pe La pelaha nt Met Ca Sea a RO hye wet ay OE a8 93 $a 8 dean’. a ance <3 1 Ninn ne Fyre 
mea i “ue ee eee a ee rd us Vege | Or es re > Tae teh hd Ah th tet aS ee Atha! Hosea a ¢ase. pai elarh ¢ fe a rN ah ates A eee ae 8 intante fame oo s8 
Sais se re i oe os a on ew ie eee ghey bred aerate Bia gray "Stacie ei late's tate cerd: bude [palareeis petites sesbeb a utinea ate adler sea eebens 
: Gate a ae We rei ee 1% ‘ ee oe eae gf fy Piast 2 a ar Sewartes a: ta apf atigd pan pets wines tap ee hens ea oO 6, Gc Ren titus ness Ac ars malin. te Soins 
bs bs ¥,' , ' bs % Ae i - 
a i g4 * 4 mm» mere Pastas d tinted ba ao vdaa Be Rois SB ena Pesee tense De ots ait: ha ie Os & be, "Oy de be SOD Mee Sa ACR eg 
Da ee A ih ont Pars 1, at de find Bey erat ep) ets tha What ta bs ER Rites whnfity Serr ert Pm bite gs be oe ne beds Gu Mageevach eh treme eh 
1 ' + t mn ay rarer "| LA aa} t + oe Pig } Arafe se ptaw & : edi pie rik Ass 9. dake 4, ,a a =P » 4+ 5 ny Pe eee Ue a Bie (shorty a-ak Ah, Bock Me FOAe de Re oa 
2 Saoeee : - . ‘ o's 2 PP fe tee ' . @ ptaty Wi % dy dod ie anon gh eregh GAR aeons ie hr ob ae wh nace gas pte bod bb ee the Oo Be alee Arh hy Band be ty te 
wee : ' ; Ae _ f out Viwwk 0 ! = hee et tal fe Avtar Pe ea piitckee pees as alae Re RYT ? nana us Mogae ii irume ae ra ach ee. Peart 
Spee 8d . 3 Wer Su , 8 : : A, tet he" , ‘Wee: rs y rr TL feces Awe a Abt Dei teh RSA tre Vee arats meson Mahotek.: ue ane a cape Bt ime kara ga Peete awed os ace at 
: i S ’ "gin "eel r pei 3-8 nf OU » pea ie be a Pd ahs 01g, of Reh Oe Be Oo Me A Oe MO pa bd ba WIE rer tee BES: $45" 1 2 Aameqee 0°47 > He bec er ee Sagi Mem Sa ts ok Rede” cea ee road 
. ek Ss : 5 UF; ' * re. eee ne ea ‘an oe oe Thee wee BATA etic s wi * Ss 2 eee a Bn AsSacte ert eel 4 *S i ieee ew bet he See Oh, As SEDER ree Biogas 
oD ieae CAs . a - 3 “ “hie Len he =o 4 ty Into Ae ren als ® Ot Nett re ade a a 
Lene : Feil " a I ‘ 2° ¢ A To ¥ rs aA 5k ee neta ne a ; kat i 44 : es Re Nit 72 ba ute} ae ath wre aa" eal # fetes etegimres Aas ay = are Atipa 
fy oeeats ' ‘ . eet na 1 § . 9 ted wm, of rae > *h eee Leeman ap ¢ : dep yet Hi Nite BG NEre a> Rhy ely oe biped eel eA mpd nhc alana v9 
y . ae a ar 2 rs "4.4 * ma Ms oe oe by * +8: ' sata 4 “ne em ee Ome mad ves sah ne yr «) eee an aa Relay At &y Sow 
ass a ae be ’ ; ve > 4 ‘ is ee Kota fedatey: baits Rete) yp ns ate vas F Rade | 4 td ary ots ns mashes ‘ Jao eg mmo ted 
. ' ‘ , ei 2 : > : K " body 28 a BS Dae ; s ’ : - i era bat 
a ' 5 eee u es ' hit oe oF Sree - ‘ hrs Wart ae Ae: Wyle athe! att tubby War®§ Preity RA Gogg de Lets awe he io ieee’ 
CeCe Cy a as Ages % kh * ore we = Ooh Le Re eas Fs ee < | Ye oted 34 54 © ee Ue a Ne ee ek ie eS Rr A phere! ore Beer 
' Ye ; pela ' tb i? i Se rte Behig tN OTN HAL Ble Samy De Ra me Me aR ate SR gy bi eee Waa . SIN A De ONT, by Apend twat Sansa) tasers 
ot Lae “A et rT ne nee . L ak = mith deur & Aa! bet eek & th tee ate os vhs am, Aus hoabbentncd Es pti at EAT Uw ON Ung Ad of $705 92% & BAD) Soe ects BA. susie aes 
. Lee ele a - . 1,¥ “ ‘ papi ew gtr , -~ madd fy , Lore es ware On 7 Ceatihe L 
id - x . oe ’ : gehen ta ’ re ls oe Se ~ hm 4! ae Sr AST * etfs arin siaresd 3e ie ays teeter tala ie Miata! Mileaarared ae Glyn peer. “—e ere 
; . Chee td . 4 * ; Fi ' . 4. > ne = ¥) Peg ahy = “i us ay ad a5 mate 4 2h What? oo Dear b dee ee ett Cer aa Y Y \: ar ids pore penaee Fast A Nate: 
; thos ‘ ’ i teas ee FH aide BD f tye em laliad } = ? ay My Fp pail wee VErTe SS asi tphad’ + Tote ms Ae Oe : Bh + Si Es m dyads is 
Z . ‘i A F - oes me day Me m ote - > na 
1 W " | a $ ay s " A i a, {19 ta) heh, eh! Ae * Sawn s Wusw ee tt wag ee a eS uate ae - wie! SCN, » lelae eines are hd eh 
aor . $ yo4 “ ae Ki 4 ie teas Cig e” Be 4 7 mie" * taran 4; tite te aor Care Te hee ret NG st ay eer = rege vl 
Ce , love soy , * ‘ ¢ ‘ ath oe “"h .slerte ge 4 By uy eae ai 4 iat ™ i: Heal pi »e) Re AP Ast 0am, se pt idesy yp & saa bey Spe Pee oh 
; ero \ ' Ae , oS: eau pa Nel Ane fya Re iy" } aye ae abe piataa?s an, Viakadace? 2 Ava LNs warpiecotice ip teyon ~T 
: a s - . : ; . Oh J i 2 a, : . he i bs 
' , er, : : : , | a ; ou more ss f aye, ‘ah OS nace ©, }! a. i A eek rll sorte orale Pets ete ~ der 373 . bak 
: ee ‘ ce a ‘ “hy 4 7h “ Mo ope Sore t. at at ea HOR as ty ae 7h Magoianetinaty ecwlargn Sarasa hone SUM as a MSC I ME are mane a Seale bnes 
: ‘ gs J ‘ : 1 tha as See. sey l# Smee os tare Oe eae 4 8D eile ie aes AREA retrial dodets rT iigaogim MAS ah Laie © bed Ad Alas aaks tos ilend beh 5 
oan : ’ Sey a ig Faritaaay fag’ ie ty Ae: & REL NE ee OB ee D Rake) th A igh oh ee <a) Deg tr Os fa ge pe Vase Oy hy Some be Ss Mette DET AE a. 
' : ‘, Od, : ae a 14), Were’ f a af # ait Ss i eo es a Bebe Sn! 2 ot oe reed ie reg P oe 
' . = : aes ; . “sg! oa i) 7a oe i seo Se*simn orcawe fue 
| 4 : ‘ rat ; "3 r $e gi &, &y 
, . ' ee oe Me a) . ; x in 8 oe as 4 on 3 hs ! wine: Aa ry saree ny Lares dase of 
- . ae @ 6 . — Bet) ie a reed - 
; DC a. Tar, A. pe ! me yu 38 ls efit aah Oy Cigar lts pe onlay oh ay 
Of ae oe ' mat og fe ys) “_ _ ety ery Hela he Hah Aecta BES Fe Bracing yess 
1 > ‘oy? > otet Aas y" or ; i . yan a pi af ata’ t's “§ 2 yt cae ithe bees pd we 
a t aan ' ? vas hae $-) ae $ ' PY ay 4 hppaaming cw Talh e 4) & 
: ; Hen vis ainmavae 4 are | TH ager as ‘eh Fei ay te een ~ » 
. : tee ~ f » 7 oy! ; of, Bh > on Ze ye i 
—— ro ae *“s" ee a aes Tee a 2A Dero ok 
ees ae : hai ; eee Ay at ee See fe kta nd al oe BF tee oF PLS 
es enue oe ar an wants 5 ee m, "Hye Se OB UM? HHeRhM 1h fies (Pas easter el ’ 
oe , ' Pate e te. on : ' = anes ' 2 ee 1% tes ; ee, A Foe . ys ie! i 
i . : ie poe i , ’ bs : ae . Aba’, Beals ge 2 webs 
, 4 _ : £ ast 7 Pe oa a ory . “4 oe cy 
. $ s soa? 4 W ey ce aby ales ect e 1h, ° r. 
% . . "by ; "a. ey" ' u a ® & He F. ’ ct 
' " cee ‘ Pe ite! te Sie Y aay ¥s 
us 13 A oe we ee ie | faa *hy'y° 4 » * eens 5 a 
of ' , qo FeM 2 : s i 4's ' * : SC i : :¥ 
' ', : . ae : " _ 14+ ihe } ye 
2 ‘ ’ : § ay *. L o fa 7 ee ae Jad slag beer ata rb. Sia: ; . 
eon : “teh 7 : 2 sot Ae the 2 sb PORN EEE ray Ries = hs 
is ' z aa" poe aes A . : a, : i, a’ te, re ae tz Poin mty Pyare) +a! Date pe a 
a we t,s - wt ‘, o fe ’ re ges Ayu £5 sh aSbls of as bE moet as af Mat 8g KB Me a RR a ae, 
cae 2 1. . ; ’ oc a r> toe : BEES ean of or attaiefetare astisia. &F rae 
aca ae , F 2 * ‘ A : Ny Ope SOT, ASE S Poa ee BM ae Sa tatie sr 2 ig 
A 2 > i } gould ‘ H og ‘OO aby ey oletala Pa Poe er OR 
2 ‘ F | ‘ ' ' se Vie4 1" ? Ya? ey Fo Suge} » id } 
: fu tune ois aut "yy? a. é A fe Tate bie te? af sive" EEN Se voy Mae hs oh 
. . ' ‘ . ' Dees et 4g i - a7? ahd + OR WEY ta anya: eats ree 3a pase oko eat 
eae a, oe i oy “a” Uh efeg fan 2 Nis ‘ 
Ho : rer: cr. Wes a. ’ ee pits Sees cs = Me rat arte a ye, Rate AMA: ; Ste, he nat 2": ok i 
fi . , ' aly : to 24 cs a rat " to J x iy alam Ga utes! Sulhl pias arith “ rie ay) hh 
. Apo: » * a an $047 sa pata’ “Byte ee oe ae * ts serie al ne sf > PR af sO ass ttl + ; a, F a 
‘ ' adv . r oe '. . % 
Fi "5 ‘ i Ae tt A | ce ' ee + "ys : tits a ie " ¢ ae : He) : TY ere ap . ao! ih 
: ‘ u 4 . r 4, ve are ae ‘ * 8 ry Psa PRN RR ‘ tal I ae E 
ae ee PE ag eg ee + 2 OR he ahh Sipe ME Ee od stata 08 Nea MI 
=i , a . ‘ oe Me Ae uae Lh wae a hy Hi , See Pen nat 3 Uy be " 1°} an Inn y abe eeont 2a ‘ Biter e yt bees we bow 
. 3 te eer o 48 ‘ La ea i ee, es gee A res Md gee ea Fle oe ES PRE age ‘ Rie’ ae"9 Ee ag ord. © 
ae one Fh reef YF) Fs } ity, Sls cra Eo A ee E 
3 pan 9 ' ‘ R tet wreak | +, . faiathogt” ol Mite ea sf 1FGS Ld ce . wa yeas ss rea ate! hog. : i ¢ 
“a an | t : ' ' Py I a , 8 2.4 ez y wpe se HR, Pia Bes ite 3 as be Sader es ay ¥ e . Pe ; rth u Nori Hy naled i AY i: fi ie “ 
: . " , ‘ : ae ; ; re 0h vere vg ' “e6 SUbatef alasas ot * fotki ote ef : Oe ee ee red ra) i wee 
28 > : ‘ ' * ? r 2 we Te tts tat 3 # s Ls PeeeY, ry » J ' pte tt tet ye s > fer: a a of ©: 9 
i Tis Fi ’ wees . Ms | - ors . er eee ee Pe rey —e a! Pa’ 5 Pe ee at Hy gee ae Be es ” or as aft! < 
" ' LY me ase ne 5 &5,% a. *, Sees eae ' » Jes Aa Saba ahaa nt it: syle ¥ g vty otf f 2) 
' : . : : 3 Ath ao OS ; - , f $7 ' #y> "fa an he 1} A os i * i ‘Y a ote pie ra Xt peice ital z tet ata fs Hi ; ye re Ear? Hees 
‘ - a: o © Es ave & e . pat? ee jaa A POMLEW Ler i 4 maf AAPL ‘vine ed A he ee 4 + 1? ote wheat at i ois at als fd 
H hi ? Hytgt 8 oye F: es ; ma ARE a t ne \° ’ vtats ra i pat 8 ag A f cia 73, € 5 tart. ry fy) suf st ot 
“ * rs . ie i tae yy ‘ . ’ 4 ort 
' 1 ‘ t's a ; ry = HD er 1} ‘ ee ‘ ' atet ny oP cee =9t ce ates rt 
; : : eigen 
Lee | a 4 . 
2 7 2 
a 8.8. . 


a> 
at 














ie 
zt 
~ 





































































































































































































' 18 
Sits,’ 
aa yop ve 
eed Sena 
. Ue bd 
'‘ “G J oa 
eet ¢ 
4 } or a '3 
. “sf 2 ih 
fe < 
— + Fy 
' 
a 
. U sigf 
oes 
: sth tee' i 
ae ¢ ey 
i Dew oft 
ai : fise 5 
; “aes shae 
: gt u ey. aes as “4 a ane i rae ac : ee whi . ne ‘ , 4 he 1x = 
: a aa : a ' : A 
y 4 3 1 ni ‘ r 
ts “4 ry ' z ; A cs a ? ' ate , 5 a “ 
fa wt , ry an ‘ ae aay ' A ‘ 7 é ° of ak 
4 ; , a? ' y i nee + ge Pite o va tie 
‘ of . o .4 ‘ s ’ e.8 ry 
oS ; 5 : ' ; r a. a 08 gee Pi 5 
— Ps Ce ae aE wise ata aig fal aah ts 
. awe ae = PR hc at 4 so 
; :. é ‘. o ve on ee i : ie," : ied o 3 ate > Bh t ‘ u 4 > Peas BET: ot eee at “ ne 
’ ? iol as es one 4, (ee ya Lk i yet Fy Sate ae a0 ai af rR ig fae pas STARE eg 
. ; a : ae a : ao Tae fz 3, n't 8 re Lp re) eae a ’ ee sak - wh ee oo ene wees 3 Se Soh 
’ tao ‘ : Oe Ue OS re ea 3 five ct Ofte 110 i™t 28 i ae] % Kier! ¥ Uae Psted 
. ve ‘ bea : ee Ten 7 te ‘ F eh te gale’ ates i : e BARE af: ; 
7 ; te ea tte Dae plana Os . ceeeeie ae Ue ete fe ignites F; ve 
. tf - yi t 4 ave ' ’ } % A , ‘ he a fi, ort dey of if}e oe Perera 
: i ey a Je fay te ie ¢ aed arts ee ae ‘ fon ’ ee Ponty wry : : : meget Hand ie 
3 . ue . . b oe iy" Pee ! . a rt , Leste . if “4 , yeeus tates . Santis os 
ae oS Uh td Rot ' . ‘aa i ae t ee ae aes ' vat oo 
wo ' s ";' t fig Pe vet a te ; a t P ee ir . 
. ' . . . . 1 ¢ ity rn Be uF 
Me 7 or rf ‘ aoe fy A eit rm vie “ie ‘a Dri 8h ee juste pp al 
‘ as is : ie le ; "es iin Mien TA Gee 
i "7 ; iG nit . i a@tr di. ae: s aaa SHR 4 ’ mer. i 
i i c ; t 3 Piss fl rifts Le es Ged ee Ie on rage Kips! Ce ee 
Ss : PY reas A hel atelite ev ffs ' nes 
Ae mi tne ? a) Cia 0a > ol ee cet ey oh ete o & Sates ates 
. tua oo » $e 2 on co ' Oe De. d Pe mon a ai Persp ott ni ° i A Spite aes 
nets la . (tee oat wm, tlie) t el Yoo Ware A, Seah ek ta Ara ee re 
: ’ i : ros ee 2 ae Fuge V1 a. opedca e eats har SEH oT oo 
. ' -. a F aa ' qe fees r ye meres fy ia Sere ste 
‘ “e : ME Let) a pt BR ea " a The pae= Pee ee Det ‘a4 
Y : a Lee G PUP O geo as Tacaea, oe LP) £ eheb itor, dt ad hana ces ete 
, ce : : Sg | OS fata i rae Gh AES as Th oa idk Frater bree ee an) i? 
. ; ‘ e = By r 3 re i f he we a FY er it : ay ¥ we rt ‘4 ie ieee Saye : 
; ' - * ’ oI Se Pie vei oval eed ld 
me , ’ $ iF ; oP ‘ ry : te i aad des Uf ry, avd iA eden 
; , os | pe Re EA Bees a ry He Gey avy posi ge aby 
; } at fli f e 7 ' ¥ i AA Fg ae Cat) a sae coed =) Maat 
, : eee os vie gs aa 4 ; fi feat “he eee ar indus 
‘ v U La 0 } Caer I UR ie peter LOU o, oCi D otw tes gf ctteae f a | cetyl tS 
ae , . ae Hire i eo hey SR to bp here vane spe ate i iii he ene. Hdd G06 
aes. Peet ah 1 ES eos ng: ritet ie Leta 0 ee pet eacaee 
= 4 "os ro - eee 2 ry) 8%, 
t 1 ‘ a 4 5 p i ° ae ‘ i Ay id ‘ fa f ates te oily net 
; ee ; a ‘en ' “ cal re ho ei cetnest rey OD ae og a ‘ tel trp ae pce : writs eRe eae Reset 
x s ’ ; re s % oF i rd x ‘ i. . 
’ uu , as ers a ' +7! Ny pap tae oy Peed Pea ay es Alot hee verte! Rec aeeon ; a ssleivie ee ree 
A Pa . Re % i ae ase FP 8 oe Fe ¢ aah Vane He PRLOY At RPRveses AF : € wagons aeasenty sas 
fs ee 4 a aR Ate ORIEN ee ites ia aerate 
: ’ ‘ af - ’ a . . RY z 
fhe j er 4 ; ’ ‘7 iar Magee 0d OOF? Fr tele ye ods peat phi 2h ile aso eeeiph 
+ 6 ‘ vr | . : . s/f or wgitnw 7 Pe? #7" Seni nentsts Pvoy Wer see “Hiker Weibtpenreew ert’ 
| . . . Dyk Peg Pea QR RE EGE Gils hit feel tdt tie tts Yiee Beite ett end edltebeeitedt taaeer Ry ae tate lagebie amen 
° i r i 9% ee PRUE wid | A, at ; * PA eg °F jist Vie paehacas veo ® wadeg ee YT ea: 2) a vadusepest gs Fie peas bea piqued) 
. t re | . ¢f L Pots cul, ty Tetaeatateny ra & eof Dopey ye wd eg ah age get ty ee A seated yy : 
. Te J " oo a8 Re aaltet rat ys He § ye Gah Glergrahtlarie Yee ti at ee 5 dyad eaee eee wie 
= 14 dq wiles FSG RG Nyd Wh he GaNNb TED Ha Meth a Metra Hkh ADR bende walnad gee ui 
: , Pu ie dias aan a UOT Les tate "WEP Saas bye ge 74, oe ate 
/ i é ( r BY vine la: k fh se oft a : i HEWES e0 i Vie ee sae than Ex 1) drags adel ge reeee pe 
i wp ring nie yi ee 8 eth & piAT ¥ ae mn thin pans ita nh Ps £0794 ge Sav: 48 geeeiegt 
jo oa¢ F i im 7% 1 * 4 iy 6 Pt ine ¥ ae ares SYR ght ai ye whi 2Y0ET GAN a Tere @isuaiet tiene 
} a 4 Vee i ee : f \ ' } b > ¥ are rear ng \3 ats $) orse Pha fh Tike 
fer t t ; j uri a at al 2M) eek add ute Raee 
ee fers Fevcies hesie! MEA Ge sanvergige Res 
: * f by Fey at ale | ge Te OE esate! t iilece ya Gita 2% 4ugs gene te 
: a : : ' LRP eRe cm AIT T Mee! te ey wpe pad hh Say Loe 
; * @.8 the ee aete Py 1 Ae Ri fae Pitta Vis hh bea tla ae erate ‘ « 
ie | ie Pept het Be the es Vat Mircette Hey cent Parana 
s ' 2? i o* Cay ie Mee Pay) Vee rae per O91 me eatery apy ts 
ji : % ed Peaitgsiemuntoaes Roh ka sk bt oi 
: site ates) Ge A , fe 145 Puy UCHR ES wee ge rors 
: r * ee 7 Vl B90 pe dup- gueyy pe yawegewED! iret 
' . st , yee OTA Nyt Br Pine ges Crgegin ery. 
‘ io 8¢ OF Fe 4 yoW PEN oS a kde Dl Guar isigyd 
: es HAM RGUNE Mein h gy Oasnacisrerg revises 
. 5 ° 134 z ev _ 1 “H G4 0 Orbs , ge Gg wre eew ange 48 xa 
1 a en one ts F yer ener eve wt fata Ket Hy nan bed space ponte bee, 
i : ny ee te Hear meted did Sa rerun estat paeree narra 
.at st ia ny 3. tabby Rares S ar si roi 4 Toeras a yD. CO atte ae RE CER D mcerEeaT Fee 
’ . is ; Pt om 4) dp Ny are Pap pe pe, ; ULC Fr i Pa i 4 be : aera pp peas PD gn 
: be ' ’ Ca cs ra IP eh ee thet Pak Ah tests oh UY BOA pe 8S aka cote) ¢ HUE ged Fe namaetesve 4 MALS Rite 
© 4 . 23 ¥ 2 si, oo ? LIS Mae hges i, 1: Qe RE ipl F PEPER. FP Ord F895 ES 9 PAUP ORT ET QS 4 se VEY RAD pi it 
’ yt 8 TINE b treg ie COM Ate PS PT PANTS Wye ae! 92 Fae feta es Pansat Sve domera iat 
7 * " aie hit Rr h err bee radon 90") ea wey eee eee 99 OTE RY GoGe Gi 
' . . . f 3 q 1? e of is PPaKACeWe ites pypiry re gi4 " TR. TAGE Neen PPD fle 





ea 














NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


SURFACE RECONSTRUCTION FROM PLANAR CONTOURS 


by 
Allan R. Jones 


June 1986 


Tnesas Advisor: Michael J. Zyda 


ApROvcamnonr PUbne rellease; distribution is unlimited 





1234913 





; 
IRITY CLASSIFICATION OF THIS PAGE 


REPORT DOCUMENTATION PAGE 
REPORT SECURITY CLASSIFICATION 1b. RESTRICTIVE MARKINGS 


JCLASSIFIED 
SECURITY CLASSIFICATION AUTHORITY 






3. DISTRIBUTION/ AVAILABILITY OF REPORT 
Approved for public release; 
distribution is unlimited 










'DECLASSIFICATION / DOWNGRADING SCHEDULE 






ERFORMING ORGANIZATION REPORT NUMBER(S) 5. MONITORING ORGANIZATION REPORT NUMBER(S) 





7a. NAME OF MONITORING ORGANIZATION 
Naval Postgraduate School 


6b. OFFICE SYMBOL 





NAME OF PERFORMING ORGANIZATION 


aval Postgraduate School ah 
: 


‘ADDRESS (City, State, and ZIP Code) 7b. ADDRESS (City, State, and ZIP Code) 


enterey, CA 93943-5000 Monterey, CA 93943-5000 









8b. OFFICE SYMBOL 
(if applicable) 






NAME OF FUNDING/ SPONSORING 
ORGANIZATION 


9. PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 








ADORESS (City, State, and ZIP Code) 10 SOURCE OF FUNDING NUMBERS 


PROGRAM PROJECT TASK WORK UNIT 
ELEMENT NO. NO. NO. ACCESSION NO. 
TITLE (Include Security Classification) UNCLASSIFIED 


Surface Reconstruction From Planar Contours 





PERSONAL AUTHOR(S) 
llan R. Jones 


TYPE OF REPORT 13b. TIME COVERED 14. DATE OF REPORT (Year, Month, Day) [15 PAGE COUNT 
aster's Thesis FROM TO 20 June 86 58 


SUPPLEMENTARY NOTATION 


COSAT! CODES 18. SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 
FIELD $U8-GROUP iietecmOonstrucmon, | Surtace Triangulation, 


a || Planar Contours 
—— |. 


ABSTRACT (Continue on reverse if necessary and identify by block number) 

iny applications of computer graphics involve the representation of a 
ree dimensional solid reconstructed from a sequence of two-dimensional 
-anar contours. Surface construction algorithms accomplish this by 
Peameeindividual pairs of contours, forming triangular surface patches, 
lich approximate the original three-dimensional solid. In this paper, 

> present an expanded algorithm which not only handles the mappings of 
ltiple contours per plane and partial contour mappings, but also 

-lows human interaction to resolve mapping problems. We include a 
-scussion’ of our algorithm's limitations and proposed solutions. 









DISTRIBUTION / AVAILABILITY OF ABSTRACT 
CXUNCLASSIFIED/UNLIMITEO (C) Same as RPT. (CJoric users | UNCLASSIFIED 

. NAME OF RESPONSIBLE INDIVIDUAL 22b. TELEPHONE (include Area Code) 
Be Michael ide | 

FORM 1473, 84MAR 83 APR edition may be used until exhausted. SECURITY CLASSIFICATION OF THIS PAGE 


All other editions are obsolete. 
UNGEASSIE LED 
EEE’ 5 2 


Approved for public release, distribution unlimited 


Surface Reconstruction From Planar Contours 
by 
Allan R. Jones 
Lieutenant Commander, United States Navy 


B.S., United States Naval Academy, 1974 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN COMPUTER SCIENCE 
from the 
NAVAL POSTGRADUATE SCHOOL 


June 1986 


ABSTRACT 


Many applications of computer graphics involve the representation of a 
three-dimensional solid reconstructed from a sequence of two-dimensional planar 
contours. Surface construction algorithms accomplish this by mapping individual 
pairs of contours, forming triangular surface patches, which approximate the 
original three-dimenisonal solid. In this paper, we present an expanded algorithm 
which not only handles the mappings of multiple contours per plane and partial 
contour mappings, but also allows human interaction to resolve mapping 
problems. We include a discussion of our algorithm’s limitations and proposed 


solutions. 
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I. INTRODUCTION 


Many applications of computer graphics involve’ nia representation of a 
three-dimensional solid reconstructed from a sequence of two-dimensional planar 
contours. These contours are obtained by some electronic sensor that records 
data from the original three-dimensional object along a finite number of parallel 
planes. The intersection between these two-dimensional parallel planes and the 
three-dimensional object forms these contours which lie along the solid’s exterior 
and interior surfaces. These contours appear as line segments on the parallel 
planes and are represented as either closed loops, open segments, or single points. 
The main purpose of surface construction algorithms is the formation of surface 
patches between these contours on adjacent planes in order to approximate the 
original three-dimensional solid. 

The problem of surface construction from two-dimensional parallel planes is 
characterized by mapping and triangulating pairs of planar contours into surface 
patches that form a display. The surface construction algorithm identifies the 
appropiate contours, including ie specific portions of those contours, that should 
be mapped. Then connections are formed by building triangular tiles between 
individual line segments from one contour and a single point from the end of a 
line segment on the other mapped contour. This tiling operation is executed for 
all the line segments in the identified contours. 


7 


Notationally, this problem has been specified as follows: 


An unknown three dimensional solid is intersected by a finite number 
of specified parallel planes. 


The only information about the solid consists of the intersections of its 
surface with the planes. Each of these intersections is assumed to be a 
simple closed curve. These curves are not completely specified; instead, a 
finite sequence of points encountered during a positive (counterclockwise) 
traversal of each of the original curves is given. The curve segment 
between two consecutive points is approximated by a linear segment, 
called a contour segment. 


We reduce the problem of constructing such an approximating surface 
to one of constructing a sequence of partial aj »roximations, each of them 
connecting two contours lying on consecutive pianes. (Figure 1.1) 


Let one contour be defined by the sequence of m distinct contour points 
PO, Pl, ..., P(m-1), and let the other contour be defined by the sequence 
of n distinct contour points QO, Q1, .... Q(n-1). We note that PO follows 
P(m-1) and that QO follows Q(n-1), and so indicies of P are modulo m 
and indicies of Q are modulo n. We wish to create a surface between the 
contours P and Q. The surface is constructed of triangular tiles between 
these two contours. The verticies of these tiles are contour points, with 
the verticies of each tile taken two from one sequence and one from the 
other. Thus, each tile is defined by a set of three distinct elements either 
of the form {Pi,Pk,Qj} or {Qi,Qk,Pj}. (Figure 1.2) 


Each tile’s boundary will consist of a single contour segment and two 


spans, each connecting an end of the contour segment with a common 
point on the other contour. [FUCHS,1977] 


This specification is mutually described in all the public literature on surface 
construction [GANAPATHY,1982] [CHRISTIANSEN,1978] [SHANTZ,1981| 
[HOGAN,1985] [FUCHS,1977]. Each of the papers uses the notation to expand 


upon the initial algorithm originally proposed by Fuchs. 


The initial action of this paper is a brief review of all previous surface 
construction algorithms, concentrating mainly on their capabilities and 
limitations. The main part focuses on the algorithm presented by Hogan. This 
algorithm is more comprehensive in that it can handle Sande contours per plane 
and partial contour mappings. ee it also does not provide a complete 
solution to the surface construction problem. Following the discussion of the 
Hogan algorithm we present a further expanded algorithm which attempts to 


resolve each of that algorithm’‘s limitations. 





Fig 1.1 - Two contours on adjacent, parallel planes. 


* - surface patch defined 


by {Qj,Qk, Pk} 





Fig 1.2 - Mapped connections into triangulated 
surface patches. 
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I]. LITERATURE REVIEW 


The method for finding an approximation by eeeealation of a surface 
defined by a set of contour lines has been the subject of past articles written by 
[KEPPEL,1975], [FUCHS,1977], [CHRISTIANSEN,1978], [SHANTZ,1981], and 
[GANAPATHY,1982]._ Each author has addressed different aspects of the 
problem. However to date, no reliable algorithm has been published which can 
successfully handle triangulating complex surfaces in all cases. The reason for this 
is that insufficient information is obtained from the contour lines regarding the 
gradients associated with the surface they describe [KEPPEL,1975]. Contour 
lines of an irregular surface, such as found in nature, do not lend themselves to 
curve fitting, or other attempts at precise mathematical descriptions 
(CHRISTIANSEN, 1978]. 

Our surface construction algorithm is based on the efforts of Fuchs, 
Christiansen, and Hogan. In order to fully understand the underlining problem of 
surface construction, a brief summary of all previous literature is presented. This 
summary focuses mainly on each of the algorithm’s capabilities and limitations. 

Fuchs algorithm for surface construction was presented in |FUCHS,1977]. His 
mobilen statement, stated in our introduction, is the basis of all subsequent 


literature. The main contributions of that paper are the concise statement of the 
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surface construction problem and a method for connecting simple. closed contonrs 
(Figure2.1 and [HOGAN.1985]). 

Fuchs’ algorithm contains three major limitations in dealing with complex 
surfaces. The first limitation is that his algorithm can only handle cases of 
simple, closed contours, with only one contour on each a the mapped planes. It 
cannot handle the more complex reer of multiple contours on adajcent planes, 
partial contour mappings. or open (non-closed) contours (Figure2.2 and 
[HOGAN.1985]). The problem with multiple contours on adajcent planes, arises 
from the fact that Fuchs’ algorithm does not provide the mechanics necessary to 
identify which of the contours should be mapped. The more general case for 
surface construction is to have multiple contours on each plane. The second 
limitation of Fuch’s algorithm is that it performs a complete contour-to-contour 
triangulation between adajcent contours. even in cases where a partial mapping is 
nore appropriate. Partial triangulation of contours 1s most often representative 
of situations in which we have dissimiliarly sized contours. The third limitation 
in Fuch’s algorithm is in its inability to handle open contours. This is the direct 
result of his algorithm’s lack of generality. A method designed to handle the 
partial contour mappings 1s also capable of handling open contours. 

In [CHRISTIANSEN,1978] we see an algorithm that is similiar to Fuch’s. 
The .major difference is a mechanism which allows human interaction to resolve 
Inapping ambiguities. This mechanism allows the user to determine the relative 
connection points in the mapping process for highly convoluted contour cases 


IHOGAN,1985].. This procedure can be quite time consuming. depending on 
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Fig. 2.1 - Triangulated pair of simple, closed contours. 





Fig. 2.2 - Example of multiple contours per plane. 
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the complexity of the data base [CHRISTIANSEN,1978]. Christiansen’s 
algorithm has the capability of handling some simple branching. Branching 
normally results from a pair of contours in one plane being mapped to a single 
contour on an adjacent plane, (Figure 2.3). This Granchine capability allows the 
algorithm to handle simple cases of multiple contours on adjacent planes. 
Christiansen accomplishes the branching capability by utilizing the following 
procedure. 
1. Introduce a new node midway between the closest nodes on the branches. 
The Z coordinate of this node is the average of the Z coordinates of 


the two contour levels (planes) involved. 


2. Renumber the nodes of the branches and the new nodes such that they 
can be considered as being one loop. (Figure2.4) 


3. Triangulate as usual. [CHRISTIANSEN,1978] 


In general this algorithm introduces a new node between the two, planar 
contours. This new node is used to form single, connected regions which are then 
processed by the original surface construction algorithm. 

The problems with Christiansen’s algorithm are its inability to handle open 
contours and its inability to Sear complex cases of multiple contours on 
adjacent planes, except through the use of expensive human interaction. 
Christiansen interestingly shuns the optimality seen by Fuchs as important by 
utilizing the heuristic of choosing the "shortest diagonal" in forming triangular 


tiles instead of minimal triangular area. As stated in his article, this heuristic is 
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Fig. 2.3 - Simple case of branching. 


tv Wan 


= ene a= _S 


introduced node 





Fig. 2.4 - Triangulation scheme for branching. 


easily implemented, fast, and works well as long as the two contours/loops are 
mutually centered and are _ reasonably similar in shape and _= 6 size 
[CHRISTIANSEN,1978]. The process considers the next two nodes of each 
contour as candidates for triangulation. After BeCerntte the lengths of all 
possible diagonals for the sarfacepaten nodal selection for triangulation results 
from the surface patch exhibiting the shortest diagonal. 

The algorithm proposed in [SHANTZ,1981] is basically an extension of Fuch’s 
and Christiansen’s algorithms. This extension includes the capability to handle 
contour defined objects which are highly branched and have holes. Handling of 


multiple contours on adjacent planes is achieved by: 


For branching contours where n contours in section i are connected to m 
contours in section 1+1, the surfaces are mapped by first concatenating the 
section 1 contours into a single large contour using a minimum number of 
minimum distance links, similarly concatenating the section i+1 contours, 
then performing the one to one mapping between the resulting composite 
contours |SHANTZ,1981]. 


Once the concatenation process is completed, Shantz uses Fuch’s closed 
contour mechanism to formulate the connections between the composite contours. 
After the connections have been formed, any extraneous connections resulting 
from the concatenation process are removed. The resolution of ambiguities 
areine from multiple contour cases requires human interaction and simular to the 


Christiansen algorithm, Shantz states that this is extremely labor intensive. 
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Shantz cites a specific case in which a set of contours from the Livingston brain 


database required many hours of contour splitting with an interactive cursor. 

The main limitations of the Shantz algorithm are its inability to handle cases 
of open contours and partial contour mappings. Additionally it can only handle 
cases of multiple contours on adjacent planes when a composite contour can be 
formed. or the ambiguities can be resolved by human interaction. 

The algorithm described by Ganapathy [GANAPATHY.1982] is a further 
improvement on the Fuchs and Christiansen’s methods of handling simple, closed 
contours. This unprovement results from using a more computationally expedient 
heuristic for triangulations |HOGAN.1985]. However, Ganapathy’s algorithm 
does not include the capabilities introduced and discussed by Shantz. Instead he 
sunply assumes a complete mapping of paired contours. which is not always the 
case. 

The problein with the Ganapathy algorithm is that it represents a general 
solution for handling only simple cases of surface construction. Capabilities for 
handling multiple contour mappings, partial contour mappings, or human 
interaction are not provided and their issues are not addressed in his presentation. 

The algorithim presented by Hogan |HOGAN,1985] is more complete than its 
predecessors in that it not only handles the simple cases of contour mapping. but 
additionally provides a more comprehensive procedure for resolving the multiple 


contours per plane and partial mapping problems. The only capabilities lacking 


from the Hogan algorithm are the one for handling branching as described in the 
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Christiansen paper and the one for human interaction for the resolution of highly 


ainbiguous mappings. 

None of the above papers provides a complete solution to the problem of 
surface construction via the triangulation of contours. What is required is an 
algorithm with capabilities for handhng multiple contours per plane, partial 
contour mappings, and which ort: simple cases of branching. In addition the 
algorithrn should provide a inechanism for human interaction for the resolution of 
highly ambiguous mappings. 

The surface construction algorithm we present handles the cases of simple 
contour mappings. multiple contours per plane, partial mappings. and in addition 
provides a mechanism for human interaction to deal with cases involving highly 
ambiguous mappings. The only capability lacking from our algorithm is that for 
handling branching as described in the Christiansen paper. <A discussion of our 


algorithm follows. with a proposed solution for handling cases involving 


branching. 
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III. SURFACE CONSTRUCTION ALGORITHM 


In the preceding section, we presented a Fee ikcion of previous surface 
construction algorithms. Here, we present a detailed discussion of our algorithm 
by first specifying the known input/output data structures. 

Surface construction of an object between a set of planar contours (Figure 
3.1) can be reduced to constructing the surface triangulations between two 
adjacent planes. The specification of the problem can be best seen by listing the 


known input data structures [HOGAN,1985]: 


* total (i) number of contours on plane i. 

* start(j,i) start of contour j on plane i. 

* length (j,i) number of coordinates in contour j on plane i. 
* type(j,i) type of contour j on plane i. 


(CLOSED LOOP, OPEN SEGMENT, or SINGLE POINT) 


* interior(j,i) value of contour j’s interior with respect to 
the contour line. 


(HIGH, LOW, or INDETERMINATE) 
* coords(XYZ,pointer,i) input coordinates for all contours on plane 
1. To isolate contour j on plane i: 


for (pointer = start(j,i) + k - 1), 
where k = 1, length(j,i). 


From the above data, we desire to produce the following output data 


structures |HOGAN,1985]): 


19 


*num coords number of coordinates generated for the two 


input planes. 


*new coords(XYZ,num coords) coordinates generated by the surface 
construction process for the two planes. 


* new conns(num coords) drawing instructions for each coordinate 


generated (SETPOINT, DRAWTO, DRAWPOINT). 


If the output data is in the form of triangular surface patches, an 


alternative data structure is required [HOGAN,1985]: 


*num_ patches number of surface patches generated for the 
input two planes. 


“new coords(XYZ) new coordinates generated by the connection 
process. 


* patches(3,num patches) a3 by num patches array of triangles. 
Our surface construction algorithm is composed of the following outlined 


steps:! 


A. INPUT AND INVENTORY COMPILATION 

The data structures defining the contours are processed to extract the 
pertinent data. This data includes the number of contours per plane, the 
coordinates defining these contours and the types of the contours. 

Additionally, two-dimensional bounding boxes are described about each 
contour for processing consideration in step 2. This compilation of data is used to 
create the data structures required for surface construction. 


' The bulk of this discussion is drawn from ,ZYDA,1984| and [HOGAN, 1985] 
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B. OVERLAP DETERMINATION AND CONTOUR ITEM MAPPING 

In this step of the algorithm, we determine which contours on adjacent 
planes have significant overlap, and which contours’ exteriors are near. This 
information is used to designate which contours ayould be connected via 
triangulations. The assignment of ouenite is accomplished through the use of a 
value for the overlap percentage. This value is computed from the areas of the 
two-dimensional bounding boxes, as seen in Figure 3.2, of each contour. The 
overlap percentage is used to give priority to contour mappings that have the 
highest percentage of total overlap area. 

In this step of the algorithm we also perform consistency checks for each 
contour pair. One such consistency check is executed using the contour interior 
specification and the overlap percentage value. Contour interior specifications are 
assigned as the value of a contour with respect to its immediate interior. As such, 
a contour is LOW valued if it is taken from the exterior of a solid object, such as 
the skin of an apple. Conversely, a contour is HIGH valued if its immediate 
interior is non-solid. Using these pieces of information, we are able to eliminate 
contour mappings of high overlap percentage which would result in an erroneous 
approximation of the original three-dimensional solid. 

To illustrate the application of this consistency check, let us consider the 
mapping example for Figure 3.3. Here we are presented with a set of contours 
taken from a solid cone standing within a hollow cone. In this case, contour 1 on 


plane 1 has a high overlap percentage with contour 2 on plane 2. 
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Fig. 3.1 - A partial set of planar contours from a 3D Z*%-orbital 
of a hydrogen molecule. 





Fig. 3.2 - Two dimensional bounding box used for determining 
overlap percentage value. 
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contour 1, plane l 









contour 2, plane l 
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contour 1, plane 2 


contour 2, plane 2 


Fig. 3.3 - Example of consistency check using item interior 


specifications with overlap percentage values. 
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However, since contour 2 on plane 2 is low valued with respect to its solid 
interior and contour 1 on plane 1 is high valued, this mapping can be eliminated. 

These values are also used to determine whether the mapping is interior to 
interior or exterior to exterior. An interior to interior mapping is one which maps 
the interior of one contour to the interior of another contour. 

This form of mapping is indicative of contours taken from a surface with a 
shallow gradient, 1. e. - a sur:ace where the mapped contours are of similar size 
and shape, and where the contours have significant overlap. An exterior to 
exterior mapping is one which maps the exterior of one contour to the exterior of 
another contour. This form of mapping is indicative of contours taken from a 
surface with a steep gradient, i. e. - a surface where mapped contours are of 
dissimilar size and shape, and where the contours overlap percentage is slight. 
Interior to interior mappings are more common. The exterior to exterior mapping 
is indicated for cases of two contours with a low percentage of overlap and 


differing interior specifications (HIGH:LOW, or vice versa). 


Ce FORM COORDINATE MAPPING FOR MAPPED CONTOUR PAIRS 
For each coordinate pair from step two, we form a complete coordinate to 
coordinate mapping. A coordinate mapping is a tentative set of triangulation 
connections between the contour pairs. There are two procedures for determining 
this initial coordinate mapping. The procedure used is dependent on the type of 


mapping found for the paired contours in the previous step (interior to interior, or 
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exterior to exterior). Additionally, both procedures try to form triangulation 
segments of shortest length, similar to the Christiansen algorithm. A general 
statement of this selection process is that we are trying to map coordinate i of 
contour n, plane 1 to coordinate j of contour m, plane 2 such that the distance 
between the two coordinates 1s ee eed An additional qualification to this 
distance minimizing criterion is that coordinate connections do not cross, i. e. - 
coordinates 3 and 4 of plane 1 are not mapped to coordinates 6 and 5 of plane 2 


respectively. 


D. CONTINUITY RECOGNITION 

The coordinate to coordinate mapping formed in step three is examined for 
continuity. Continuity, in this case, is defined as follows. First, we form 
continuous sets of coordinates from the coordinate mapping such that each 
coordinate of each set is constrained within a coordinate tolerance and within a 
distance range. The coordinate tolerance factor is a ratio of the number of 
coordinates in the larger contour divided by the number of coordinates in the 
smaller contour times a window value. The tolerance factor is used to group 
coordinates into a single set based upon their mapped coordinate number being 
within plus or minus tolerance of the last mapped coordinate added to the set. 

The tolerance sets formed are then compared for overlapping distance 


ranges. Any sets that have overlapping distance ranges are then merged. The 
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merged set with the smallest distance in it is the set of coordimates for which 


connections should be generated. All other coordinates are left unconnected. 


E: MAPPING CANCELEA TION 

Once we have decided to generate the connections for a part of a contour, 
we cancel any further mappings to that piece of the contour. This operation is 
required for partial mappings in which two or more contours on one plane are to 
be mapped to a single contour on another plane. Also, this cancellation precludes 


connecting contour points which have already been selected for connection. 


io CON SEC TOs FORMATION 

We generate the coordinates for the triangulation connections specified in 
step four. "In between" coordinates. coordinates not directly mapped but within 
the tolerance factor for the connection Mapping, are also added to the picture. 
The goal of the process is to form: minimum area triangular surface patches for 


each segment of the inapped connection region. 


CG. EDIT CONTOURS 

We extract the contour coordinates from the input data file and use them 
to create contour defined objects. These contour defined objects generated are 
then available for the user to remove or save for reprocessing by the surface 


construction algorithm as necessary. 


lal RELAX HEURISTICS 

In this procedure, we allow the user to input his own values for the three 
heuristic values (overlap percentage minimum, boundary tolerance percentage, 
tolerance multiplier) utilized by our surface construction algorithm. The user has 
the option of changing one or all three. Once these values have been entered. the 
information is used in the connection process of our algorithm to produce more a 
correct mapping between the planar contours. 

iN Input and Inventory Compilation 

The input data to the algorithm consists of the contour descriptions 

for two adjacent planes of a three-dimensional solid. The purpose of this step of 
the algorithm is to segment this data into separate contour descriptions and to 
determine the individual characteristics of each contour. Figure 3.4 consists of 
two adjacent planes, each having three concentric rings of similar shape and 
continuity. Figure 3.6 consists of two closed loops on each of its planes. Plane 1 
has two small interior lobes, while plane 2 has one large surrounding contour with 


a small interior contour. The contour descriptions for these figures are composed 


of: 


- the starting coordinate location, 
- the total number of coordinates, 
- the contour types, 


- the interior values, and 
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- the contours’ two-dimensional bounding boxes. 


With the exception of the interior values, all of these characteristics are easily 
obtainable from the input data. 

As noted, the contour interior specification is the only piece of data 
which is difficult to obtain. It reamires an evaluation of the data values lying 
along and interior to the contour (see Figure 3.3). If these values are not 
contained in the input data, a mechanism is provided to allow for user 
specification of contour interior values. The range of interior values is HIGH, 
LOW or INDETERMINATE. The problem that occurs without this value 
concerns the contour pairing problem encountered in multiple contour situations 
where contours are closely spaced and of similar shape. Here, some form of 
human interaction is necessary to designate which pairs of contours should be 
mapped together. If an interior value is not available, and the mapping situation 
is not complex, it can be set to INDETERMINATE without surface construction 
degradation. 

De Overlap Determination and Contour Mapping 

The overlap determination and contour mapping procedure of the 
surface construction algorithm is the process by which tentative contour to 
contour mapping assignments are made. The contour characteristics which are 
necessary for this procedure are the two-dimensional bounding boxes and the 
contour interior specifications. This mapping process is the key component in the 


disambiguation of multiply paired contours. 
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Fig. 3.4 - Example of multiple contours per plane on adjacent 
planes. 
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Fig. 3.5 - Connection of Figure 3.4. 
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Fig. 3.6 - Example of a set of contours requiring partial mappings 


and an exterior to exterior mapping; (1,1) and (2,1) to (2,2). 
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Fig. 3.7 - Connection of Figure 3.6, with contour interior values 


for each contour. 
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The overlap determination and contour mapping procedure is 
accomplished in the following manner. First, the two-dimensional bounding box 
of each contour on plane 1 is compared for overlap with the two-dimensional 
bounding box of each contour on plane 2. The ACO AS which define these 
bounding boxes are the minimum and maximum X and Y coordinates from each 
of the contour descriptions. (Additionally, these coordinates are adjusted by a 
constant value to promote overlap for exterior to exterior mapping situations.) 
From this operation, a table called the overlap table is produced. It is a two- 
dimensional table that contains a value for each possible pairing of contours 
between the two planes. The value recorded in each table entry indicates the 
extent to which each contour overlaps. If there is no bounding box overlap for a 
pair of contours, a value of 0.0 is recorded in the table. If there is overlap, the 
value recorded in the table represents the percentage of overlap with the larger of 
the two contours. This value is computed by dividing the area of the bounding 
box overlap by the area of the bounding box of the larger contour. 

After the overlap percentage has been computed for a contour 
pairing, it is used in conjunction with the interior specifications to determine the 
mapping type for the contour pair. An interior to interior mapping is indicated 
when a high percentage of overlap (greater than 10%) exists for a pair of 
contours. A consistency check for matching interior specifications is performed for 
every pair of contours that exhibits this high an overlap. The consistency check 


requires that each contour pair have either HIGH:HIGH, LOW:LOW, or 
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INDETERMINATE:anything (HIGH or LOW) interiors. Contour pairings with 
high overlap but inconsistent interior specifications result in an adjustment to the 
overlap table of 0.0 percentage of overlap. An exterior to exterior mapping is 
indicated when the overlap percentage is low (less than 10%) and item interiors 
are non-matching. Finally, al contenns with low overlap percentages and 
matching interiors are zeroed in the overlap table. 

Figures 3.8 and 3.9 graphically represent the overlap determination 
and contour mapping for Figures 3.4 and 3.6. Included in these figures are the 
overlap tables produced by this procedure. The table in Figure 3.8 shows three 
valid overlap percentages for three different contour pairs: (1,1) - (1,2), (2,1) - 
(2,2), and (3,1) - (3,2). Four of the entries have been zeroed by the consistency 
check mechanism. Without this capability, high valued overlap percentages 
would appear in the overlap table with human interaction required for their 
disambiguation. The table in Figure 3.9 shows two high overlap percentages and 
two low overlap percentages. This data indicates that contours (1,1) and (2,1) 
both map interior to interior with contour (1,2). The low overlap percentages 
indicate that contours (1,1) and (2,1) map exterior to exterior with contour (2,2). 

Sr Form the Coordinate Mapping: Interior to Interior 

The coordinate mapping formation procedure for each coordinate 
pair having a non-zero overlap (in the overlap table) begins with the pair having 
the largest overlap percentage. All remaining steps in the surface construction 


algorithm are carried out on this pair before the next pair of contours is 
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Fig. 3.8 - Bounding boxes and overlap 
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Fig. 3.9 - Bounding boxes and overlap table produced for Figure 3.6 
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considered for mapping. Mapping paired contours is on a largest to smallest 
overlap percentage criteria. Since exterior to exterior mappings are indicated only 
in situations where the overlap percentage is low, they are considered for mapping 
only after all interior to interior mappings have een performed. This study 
follows that ordering and aompictes the description of the interior to interior 
mapping process before considering the separate process necessary for exterior to 
exterior Mappings. 

The first operation performed on an interior to interior overlap pair 
is the determination of which contour is interior to the other. This assignment is 
accomplished by comparing bounding box areas for the contour pair and 
designating the contour as interior with the smaller area. Once the interior 
contour assignment has been made, the center coordinate of that contour’s 
bounding box is computed. 

‘The knowledge of the center coordinate of the interior contour is 
used in the following manner. For each coordinate of the inner contour, we 
determine which coordinate of the outer contour is closest to a vector drawn from 
the center coordinate of the inner contour through the coordinate of the inner 
contour (see Figure 3.10). We add the qualification that the outer coordinate 
selected by this procedure must be farther from the center coordinate than the 
inner coordinate. Also, the outer coordinate must be on the same side of the 
vector as the inner coordinate. The outer coordinates selected by this mapping 


process are recorded as the tentative coordinate map coordinate for each inner 
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Fig. 3.10 - Vector radiating from center coordinate through the 


interior coordinate towards the outer contour for tentative mapping 
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Fig. 3.11 - Example of a case where tentative mapping coordinates 


and associated distances vary greatly. 
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coordinate. We also record the two-dimensional distance from each inner 
coordinate to its tentatively mapped outer coordinate. The resulting data 
structure contains the mapped outer coordinates with the distance to the inner 
coordinate to which it is mapped. 

The tentative connection map for Figure 3.4 is very good. Due to 
the similarity in size and shape of the mapped contour pairs, there is very little 
variation in the mapped distance values and the coordinates selected for mapping 
appear sequential. On the other hand, it can be seen in Figure 3.11, that large 
variations in distance values result from this tentative mapping process, and 
mapped outer coordinates appear with large gaps in the sequencing. This is due 
to the dissimilarity of the contour pair; the inner contour is relatively simple and 
much smaller than the convoluted outer contour. The procedure used to delineate 
a correct mapping from this tentative mapping is described below. 

a. Continuity Recognition 

The continuity recognition procedure uses the tentative 
connection map and associated distances for a pair of contours to determine the 
set of coordinate mappings that should be made for that pair. In the previous 
step of the algorithm, we produced the tentative connection map for all of the 
coordinates of the inner contour. This provides a rough approximation of the 
final mapping, but it must be noted that all of the inner coordinates may not 
necessarily be involved in the final mapping for that pair. The continuity 


recognition procedure builds sets of coordinate mappings that are both continuous 
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and of similar mapped distance range. These continuity sets are then used to 
determine the coordinate sequences that should comprise the final connection 
niapping. 

The first step in this procedure is to assign each coordinate 
pairing of the tentative connection map to an initial continuity set. This is 
accomplished by stepping through the coordinates of the inner contour in 
sequence and comparing each coordinates’ mapped outer coordinate to the last 
coordinate added to the last created continuity set. If that coordinate is within a 
tolerance factor of the last coordinate added, it is added to that set. If the 
coordinate in question is not within tolerance, a new set is created with that 
coordinate mapping as its start. The tolerance factor used is a ratio of the 
number of coordinates in the outer contour divided by the number of coordinates 
in the inner contour times a window value. (The window value will be discussed 
in the next chapter.) 

To illustrate this continuity set assignment, let us refer to the 
examnple in Figure 3.11. Here. the tolerance factor is 10 and the last coordinate 
considered was inner coordinate number 24. The next coordinate considered is 
coordinate 25, which 1s mapped to outer coordinate 53. This coordinate is within 
the tolerance factor of 10 and is added to the last created continuity set. Inner 
coordinate number 26 is mapped to outer coordinate 69. 

This outer coordinate is well outside of tolerance with the 


last coordinate added and therefore, a new continuity set is created with this 


coordinate Mapping as its start, 
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This initial step of the continuity recognition process is a fast 
method for aggregating coordinate map pairs. In addition to building the initial 
continuity sets for the tentative maping. we keep track of the minimum and 
maximum mapped distances for each continuity set. These values are used for 
merging continuity sets in the next step of the process. 

The initial sets generated for Figures 3.4 and 3.6 are of 
particular iuterest. This step of the continuity procedure placed all of the 
tentative mappings for the coordinate mapping pairs for Figure 3.4 into a single 
set. This can be attributed once again to the contours) simmlar shapes aud sizes. 
On the other hand. coordinate mapping, pairs for the mapping (LE) - (1.2) of 
Figure 3.6 resulted ino initial continuity sets with varying distance ranges (Figure 
oi 2. 

Once the initial continuity sets have been created for a 
contour pairing. we merge any sets that have overlapping mapped distance 
ranges. This merging process reduces the total nurmber of sets and further 
aggregates the coordinate pair mappings to. sets with coordinate number 
continuity and distance range similarity. In reference to our examples, uo 
continuity set merge was required for Figure 3.4. due to its singular imitial 
continuity set. Figure 3.12 shows the inital sets with distance ranges and the 
merged sets with distance ranges for the contour pairing (1,1) - (1.2) of Figure 
3.6. It is shown that the 35 initial continuity sets have been merged into 3 sets of 
non-overlapping distance range. 
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Total Initial Sets Total Merged Sets = 


Min. Set Min. 
Dist. Name Dist. 
0.0176 1 0.0176 
0.1769 Ze 0.1769 
0.6067 3 0.6067 
0.1769 

0.0176 





Fig. 3.12 - Initial continuity sets and merged continuity sets for 


contour pair of Figure 3.6. 
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Fig. 3.13 - Bounding box overlap for exterior to exterior mapping. 


Only the coordinates within the overlap area are mapped. 
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After we have merged continuity sets, we need to determine 
which of those sets of coordinates mappings is the one that should be used for 
connection formation. The choice is clearly the set aan the smallest distance 
range. With this decision, we validate all coordinate pairings that are members of 
this smallest distance set, and cancel all other coordinate pairings for that set of 
contours. 

b. Mapping Cancellation 

The validated coordinate connection map for the contour 
pair has significance beyond indicating which coordinates need to have connection 
segments generated. It also indicates "filled" connection positions. By filled we 
mean that once we have formed connections to a coordinate segment of a contour, 
that segment should not be reused for any further mapping that occurs for the 
two current, adjacent planes. This mapping is both checked and recorded at this 
stage of the algorithm. Mapping cancellation examines the coordinate mappings 
for which a validated mapping has been assigned. If either of the two 
coordinates. inner or outer, has been assigned to a higher priority mapping for 
this pair of planes, then that mapping is cancelled. Once these connections have 
been struck from the connection map, all remaining validated connections are 
recorded as filled. 

An additional tasking of this cancellation process concerns 


whether the mapping of either contour resulted in all coordinates defining that 


41 


contour being included in the mapping. In that case. all other possible pairings 
with the completely mapped contour are cancelled. This is accomplished by 
zeroing the overlap on that contour’s row or column of the overlap table. 
ee Connection Formation 

When the Aor: steps have been completed for a pair of 
contours, the remaining process of generating the appropriate line segments is 
relatively simple. The final coordinate mapping for the inner contour is examined 
for continuous segments of validated connections. When a continuous segment is 
defined, the beginning and ending coordinates of that segment (for both the inner 
and outer contours) are used as boundary pointers for connection formation. The 
coordinates in between those pointers are stepped through one at a time by a 
process whose purpose is to generate the minimum area triangular surface patch, 
as defined in our introduction. The surface patch is formed by using a line 
segment from one contour as the triangle’s base, and a coordinate from the other 
contour for the triangle’s third point. The minimum area selection is 
accomplished by a procedure that chooses the next line segment between the 
contours that is both the shortest and within the mapping specified for the two 
contours. This is identical to the heuristic used by Christiansen in 
[CHRISTIANSEN,1978]. Differing coordinate rates between the two contours are 
taken care of by using the coordinate ratio (from the continuity tolerance factor) 
between the contours. This ratio allows the process to generate several line 


segments emanating from a single coordinate where there is a coordinate rate 
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differential between two mapped contours. The lines generated by this procedure 
for Figures 3.4 and 3.6 are shown in Figures 3.5 and 3.7, respectively. 
4, Form the Coordinate Mapping: Exterior to Exterior 

We begin the exterior to exterior cea at the same point 
of the algorithm where we departed in the description of the interior to interior 
mapping process. In keeping with our ordering criteria for mapping contour pairs, 
we examine the contour pair requiring an exterior to exterior mapping which has 
the highest overlap percentage in the overlap table. All remaining steps of the 
algorithm are carried out on this pair before the next pair of exterior to exterior 
contours, in largest to smallest overlap area, is considered. 

In Figure 3.13, we are presented with an enlarged view of the 
bounding box overlap area of the contour pairing (1,1) - (2,2) of Figure 3.6. This 
area of overlap contains all of the coordinates from both contours which will be 
involved in the connection mapping. The first operation performed on an exterior 
to exterior mapped overlap pair is the determination of the set of coordinates in 
both contours that is within the overlap area. The contour with the smaller 
number of coordinates in the overlap area is used in the formation of a connection 
mapping between the contour with the larger number of coordinates in the 
overlap area. The basis for this connection map is the determination for each 
coordinate (in the smaller coordinate set contour) of the coordinate in the other 
contour coordinate set that is the shortest distance away. This determination is a 


simpler version of the distance minimizing process for connection set assignment 
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of interior to interior mappings. The product of this process is the connection 
map for the pair of contours. The use of continuity sets is not necessary for 
exterior to exterior mappings due to the relatively small number of coordinates 
which comprise the connection set. 

Once we have ic this connection set, we use the same 
mapping cancellation and connection formation procedures as described for the 
interior to interior mappings. The connection formation procedure again uses the 
connection set mapping to find continuous segments of validated coordinate 
assignments. The continuous segment thus defined is used to form triangular 
surface patches for all line segments and coordinates within that segment. The 
final connection formation for the exterior to exterior mappings, (1,1) - (2,2) and 
(2,1) - (2,2) of Figure 3.6, are shown in Figure 3.7. 

5. Edit Contours 

The purpose of the edit contour process is to allow user interaction 
in identifying the planar contours that pose a problem for our surface construction 
algorithm. The contour coordinates are obtained from the input data file and 
used to create triangulated surfaces. Once the triangulated surfaces are 
generated, we can utilize the picking mechanism of the IRIS-2400 graphics system 
for editing. 

With this process, the user can remove contours that produce valid 
connections. The user can then concentrate his efforts on the contours that 


produce invalid results. After the problem contours are identified and selected by 
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the user, they are saved to a file for later reprocessing by our algorithm. The user 
can then recall the file containing the problem contours and in conjunction with 
the "relax heuristics" procedure possibily force a valid connection by adjusting the 
heuristic values used by our surface construction Oren 

6. Relax Heuristics 

The purpose of the relax heuristics procedure is to allow the user the 
option to adjust the three heuristic values used by our surface construction 
algorithm. By adjusting these values, connections between contour pairs that 
might otherwise be disregarded can be possibly coerced. 

The first heuristic value is the overlap percentage minimum. Step 
two of our algorithm determines the percentage of overlap between contours on 
adjacent planes. These percentages are used as a consistency check for matching 
interior specifications. We apply our overlap percentage heuristic in the final 
phase of this pairing procedure. Contour pairs having an overlap percentage 
minimum, with matching interior specifications, are mapped interior to interior. 
Contour pairs having non-zero percentages below the minimum, with non- 
matching interior specifications, are mapped exterior to exterior. All other 
contour pairs are disregarded. 

The value that is preset in our algorithm for the overlap percentage 
minimum is ten percent. This value, through experimentation, results in the 
greatest number of correct contour pairings. However, some contour pairs which 


should be mapped are disregarded because of this selection for the overlap 
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percentage minimum. Figure 3.14 is an example of such a situation. In that 
figure, we have a pair of contours with matching interior specifications 
(HIGH:HIGH), and in addition having an overlap percentage of less than ten 
percent. By our preset overlap percentage minimum saline, this contour pair is 
not considered for mapping and emia unconnected. But by allowing the user 
to adjust the overlap percentage minimum for an occurrence such as seen in 
Figure 3.14, an appropriate connection can be generated. 

The second heuristic value is the boundary tolerance percentage. In 
the initial two steps of a0e algorthm we determine the contour item two- 
dimensional bounding box values and then use them for overlap determination. 
Instead of creating the bounding box from the minimum and maximum X and Y 
coordinates, we adjust the bounding box values by a percentage in order to 
promote mappings. If only the minimum and maximum X and Y coordinates 
were used to describe bounding boxes. situations such as seen in Figure 3.15 
would go unconnected. In that case, we see a bounding box created from the 
minimum and maximum X and Y coordinates. This results in a zero percentage 
overlap and no connections are generated. This is an unsatisfactory situation 
since the contours should be mapped. By allowing the user to adjust this 
heuristic value, opportunities are now available for user intervention to handle 
mapping situations that would otherwise be neglected by our algorithm. 

The last heuristic value is the tolerance multiplier. When handling 


an interior to interior mapping, our algorthm utilizes a tolerance factor for the 
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Fig. 3.14 - Example of a contour pair which should be mapped, but 
would be disregarded due to overlap percentage below the minimun. 





Fig. 3.15 - Example of contours’ 2D bounding boxes created strictly 


from the min and max X and Y coordinates. Resulting overlap = 0. 
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determination of the initial continuity set assignments. This tolerance factor is 
based on a ratio of the number of coordinates in the outer contour divided by the 
number of coordinates in the inner contour times a window value. This window 
value is a constant value used for the selection a appropriate mapping 
connections. Again by the ieee user to adjust this heuristic value, we provide 
opportunities to handle mapping cases that might otherwise not be included by 
our preset value. 

We have presented a thorough discussion of our algorithm for 
surface construction. Particular attention has been devoted to the strengths of 
our algorithm, specifically its capabilities for handling multiple contours per 
plane, partial contour mappings, editing contours, and relaxing the heuristic 
values. This algorithm has proved to outperform all previous algorithms in 
surface construction via the triangulation of contours. In addition, with the 
incorporation of the edit contour procedure and the heuristic relaxation 
procedure, our algorithm can solve mapping situations that would otherwise be 
neglected. Although we have provided more capabilities for our surface 
construction algorithm, we still have some limitations. In the next chapter. we 


address those limitations. 
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IV. ALGORITHM LIMITATIONS 


In the previous chapter, we discussed the capabilities of our surface 
construction algorithm, Ppheaving its handling of multiple contours per plane 
and partial contour mappings. Additionally, we described its newest feature of 
providing user intervention in editing contours and adjusting the three heuristic 
values utilized by our algorithm. However, there still exist contour mapping 
situations which our surface construction algorithm can not handle. These 
situations together with suggested solutions are described below. 

The first mapping situation involves simple branching of one contour on 
one plane to two or more contours on an adjacent plane (Figure 2.3). When 
presented with this case, our algorithm produces an incomplete contour mapping 
because of missing data. Our suggested solution to this problem is based on a 
concept described in the Christiansen paper [CHRISTIANSEN,1978]. A 
procedure could be created to introduce a new node between the closest nodes of 
the branches. The Z coordinate of the new node would be the average of the Z 
coordinates of the two contour levels involved. Once the new node is in place, 
triangulating as usual will produce the desired contour mappings (Figure 2.4). 

The next mapping limitation occurs in situations where highly convoluted 
contours, with extreme narrowings, are mapped interior to interior. The problem 
with this mapping situation comes from our algorithm’s interior to interior 
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dependence on the overlap region bounding box’s center coordinate for the 
tenative coordinate mapping. For the section of the contour near the coordinate 
center, where the center coordinate is central, the tentative coordinate mappings 
are fairly good. However, for the section of the contous on the other side of the 
narrowing, where the center eoond inate is no longer central, the tentative 
coordinate mapping is erroneous. The limitation comes when the tentative 
mapping is so bad that the continuity recognition procedure fails. This causes the 
contours to be incorrectly left unconnected [HOGAN,1985]. 

Our solution to this situation is relatively simple and within the scope of 
our algorithm. Segmenting the convoluted contour at the extreme narrowings, 
allows treatment of each open segment of the convoluted contour as a separate 
entity. By utilizing our existing algorithm, we can produce new centers for these 
separate contours and thereby generate coordinate mappings. These mappings 
will result in a better approximation of the original object. To incorporate this 
capability into our present algorithm would only require a means for partitioning 
the convoluted contour. This partitioning method can be achieved either through 
user intervention or through some automatic mechanism. 

The next mapping limitation also deals with interior to interior contour 
Mapping situations. In cases where sections of a contour are closely parallel with 
the connection vector drawn from the center coordinate of the inner contour, 
erroneous mappings are produced (Figure 4.1). Appropriate connections are 


generated for segments of the outer contour which are nearly perpendicular to the 
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tentative connection vector; however, the tentative connections start to falter as 
the contour segment nears parallel with the connection vector. 

The same solution recommended for handling highly convoluted contours 
with extreme narrowings will correct this problem. The quality of the tentative 
coordinate mapping can be Pilar by partitioning the original contour 
into open segments and mapping them separately. 

The final limitation concerns an interior to interior mapping in which the 
inner contour is not contained in the outer contour. This situation is indicative of 
contour data taken from a toroidal object. The limitation of our algorithm in this 
case is caused by using a tentative connection vector originating from the center 
of the inner contour. Since the two contours are not mutually centered, the 
displacement between the two center coordinates results in only generating 
mappings for that section of the outer contour which is on the same side of the 
tentative connection vector (Figure 4.2). The end result is a partial mapping of 
the two contours which really should be totally connected. 

Our suggested solution to this mapping problem is again based on a 
concept described in the Christiansen paper [CHRISTIANSEN,1978]. For this 
situation, Christiansen recommends a translation procedure onto a unit square, 
centered at (0,0). The idea behind this procedure is to translate the two contours 
in such a way that they become mutually centered within the unit square. Once 
translated, our interior to interior algorithm would produce the desired tentative 


mappings for the contours’ original coordinates. This procedure would then allow 
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Fig. 4.1 - Example of situation resulting in an erroneous tentative 
coordinate mapping where contour segment becomes near parallel 
with the tentative connection vector. 
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Fig. 4.2 - Example of a situation where two contours are mapped 


interior to interior which would result in an incomplete mapping. 
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the appropriate connections to be formed in the final step of our surface 
connection algorithm. 

It has been our purpose in this chapter to discuss the limitations of our 
surface construction algorithm and provide our mareectcd solutions. We contend 
that our algorithm resolves the mites contours per plane and partial mapping 
problems. Additionally, with the added features of contour editing and heuristic 
relaxation our algorithm can handle mapping situations that would otherwise be 
neglected. However, we must concede that our algorithm is not a total solution to 


the surface construction from planar contour data problem. 


ood 


V. CONCLUSION 


It has been our purpose in this paper to present a expanded algorithm for 
the surface construction of a three-dimensional object from a set of its planar 
contours. The main thrust of this paper has been devoted to the capabilities of 
our surface construction algorithm. Specifically, our algorithm’s ability in 
handling multiple contours per plane and partial contour mapping problems as 
well as user interaction procedures for editing contours and relaxing heuristics 
have been presented. Additionally, we identified the limitations of our algorithm 
and discussed our proposed solutions for these problems. 

Although we have expanded our algorithm beyond what was presented by 
Hogan |[HOGAN,1985], we still have not provided a complete solution to the 
contour mapping problem. Further work is needed to resolve the limitations of 
our surface construction algorithm as described in Chapter IV. It is quite possible 
that the corrections of the limitations identified will not yield a complete solution 
to the contour mapping problem. However, their rectification will greatly enhance 


our algorithm’s ability in handling surface reconstruction from planar contours. 
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